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ABSTRACT 

The general Lagrangian containing the couplings of the Higgs scalars to Majorana 
neutrinos is presented in the context of singlet Majoron models with intergenerational 
mixings. The analytical expressions for the coupling of the Majoron field to fermions are 
derived within these models. Astrophysical considerations imply severe restrictions on the 
parameters of the model if the singlet Majoron model with three generations is assumed 
to be embedded in grand unified theories. Bounds that originate from analyzing possible 
charged lepton- violating decays in terrestrial experiments are also discussed. In particular, 
we find that experimental searches for muon decays by Majoron emission cannot generally 
be precluded by astrophysical requirements. 
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Astrophysical considerations play an important role in constraining the strength of 
the coupling of Nambu-Goldstone bosons to the matter [|l|. Among the various types 
of such extraordinary light particles (e.g. axions, familons, etc.) [H, which are acossiated 
with the spontaneous breakdown of some global symmetry, Majoron J° is a massless pseu- 
doscalar boson arising from the breaking of the baryon— lepton {B — L) symmetry 0-0. 
In such scenarios, apart from the Standard Model [SM) Higgs doublet $, an SU (2)®f/(l) 
singlet S is present which gives rise to AL = 2 Majorana mass terms tum when S couples 
to right-handed neutrinos 0. In general, models with right-handed neutrinos can natu- 
rally account for possible lepton- number violating decays of and particles @]-[|n[ 
induced by Majorana neutrinos at the first electroweak loop order. The non-decoupling 
physics that the heavy Majorana neutrinos can introduce ^ [11]] lead to combined con- 
straints both on neutrino masses and lepton-violating mixings. Similar non-decoupling 
effects occur when one considers Majoron couplings to fermions. On the other hand, Ma- 
joron models can naturally be embedded in grand unified theories {GUT) like the S'O(IO) 
model [|12[. In such models the B — L scale is determined by the vacuum expectation 
value {VEV) of the singlet scalar, i.e. < S >= w /\/2 and the neutrino Dirac mass matrix 
mo may be related to the u-quark mass matrix Mu by mr, = Mu/k, where k = 0{1) 
represents the running of the Yukawa couplings between the GUT and the low-energy 



scale [|T3|. There are also scenarios where m^j can be proportional to charged lepton mass 



matrix Mi [|T^, |T3|. Since the B — L scale strongly depends on the mechanism that the 
5*0(10) gauge group breaks down to f/(l)em [05 "we will treat if as a free parameter 
of the theory that should be constrained by our forthcoming considerations. The light 
neutrino masses m^^ can generally be estimated by 

1 r 

m^^ ~ - rriD m^. (1) 

rriM 

The mass hierarchical pattern mentioned above should also hold for the heaviest family 
which implies that the biggest eigenvalue of m^) will be in the range between 10 and 



100 GeV for 100 < rrit < 180 GeV [|T^. If astrophysical constraints are imposed on the 
J° — e — e, — u — u and — d — d couplings, one derives the bound on 

< $ > 



tan (3 



< S > 



V 

w 



< 10' 



(2) 



for models without interfamily mixings. The situation becomes more involved if mixings 
between families are considered. This realization will be illustrated by the present work. 



We will finally discuss the bounds from low-energy experiments [|15], on the off-diagonal 
coupling of the Majoron to two different charged leptons. 

First, let us briefly describe the low-energy structure of the singlet Majoron model. 
The scalar potential of this model which should be non-trivial under the f/(l)y group is 
given by g, |] 



+ 5(<l>"f<l>)(S^S). 



(3) 



If all the stability conditions in this model are satisfied (i.e. Ai, A2 > and A1A2 > 5) p, 
the above potential can always be minimized by the following Higgs-field configurations: 



G+ 

V ^72^ 



and 



w 



S = ^ + 



a 



V2 



(4) 



V2 J 

After the spontaneous breakdown of the SU{2)i ® U{1)y gauge group and diagonalizing 
the Higgs mass matrix, one obtains two CP-even Higgs fields (denoted by and S^) 
and one massless CP-odd scalar, the Majoron J°, while the would-be Goldstone bosons 
G+, G° give simply mass to bosons, respectively. The weak eigenstates 0° and 

(T° are related to the corresponding physical mass eigenstates through 



^ cos^ -sine \ ( H^^ 

co 



•smQ cos 6 



(5) 



where 



tan29 



25 tan /3 



A2 — Ai tan^ f3 

The Yukawa sector containing all the relevant Higgs couplings to neutrinos reads 



TOr. 7/ ^ + 77° 7/ ^ 



+ - u^R niM J^R 



1 - z J° 



(6) 



(7) 



In Eq. (7) we have assumed the absence of Higgs triplets ^7j, which seem now to be ruled 
out by the LEP data on the invisible width. The interactions of J°, and with 
the 2nG Majorana neutrinos rii - hq denotes the number of generations - are generally 
described by the following Lagrangians: 



int 



n 



■J 5^ 



H 

int 



9w 



2{ce - set/s)^ 



(9) 



'-'mt 



9W 



fm„. + rrin ) f ReCj, ^_Ji — ^ 



rij , 



(10) 



where we have used the abbreviations = sinx, = cosx, = tana; and defined 



no 



C.J = E u;:,u, 

k=l 



kj 



The 2nG x 2nG unitary matrix If^ is responsible for the diagonalization of the 2nG x 2nG 



neutrino mass matrix M'^ which is of the "see-saw" form HTH 

/ 



rriD 



(12) 



\^ mj^ niM J 



The first eigenvalues of are identified with the ordinary hght neutrinos, z/g, v^, 
etc., whereas the remaining hq Majorana neutrino states are new particles provided in this 
model and should be heavier than the boson to escape detection at LEP experiments. 
The 2nG neutral leptons rii are related to their weak eigenstates ^. and v^r. through 
the following unitary transformations (assuming the convention of summation for repeated 
indices) : 



V 



\ 



oc 



■3 ' 



V 



oc \ 



(13) 



It is now easy to see that in the limit of tan/5 and 0^0, the fields and J° 
decouple fully from matter and only one Higgs field, H^, couples to Majorana neutrinos. 
This scenario has explicitly been described in , where for our purposes we will repeate 
here the interactions of the and boson with the Majorana neutrinos. They are 
given by the Lagrangians: 



■'int 



'-'mt 



9w 



2v^ 



W ^' k 5i^j7^(l -75) + h.c. 



9w 



-Z"'^ n,7^[zIm(Q,) - 75Re(a,)] 



(14) 



(15) 



4 cos 9w 

Moreover, the couplings of the charged would-be Goldstone bosons are written down 



int 



9w 



-G li [mi^Bi^j{l - -f^) - Bi^j{l + j5)m„^] rij + h.c, (16) 



2V2MW 
where 

no 

Bia = EvlkU;:; . (ir) 

k=l 

In Eq. (17) is a unitary matrix relevant for the bi-diagonalization of the charged lepton 
mass matrix MK 

At this point it is important to mention that the presence of Majorana neutrino 
interactions in the Lagrangians (8)-(10) violates the CP symmetry of the model. In 
particular, the fact that H^, 5° and J° couple simultaneously to CP-even {'.niUj) and 



CP-odd {■.riii'j^nj) operators gives rise to CP-violating transitions between states with 



different CP-quantum numbers [20|. For example, one finds a non-zero contribution when 
computing selfenergy graphs induced by Majorana neutrinos between the CP-even Higgs 
fields and the CP-odd states Z^, J°. All these transitions turn out to be propor- 

tional to the CP-odd combinations ImC^^^- which change sign when a CP conjugation is 
applied to the vacuum polarization terms. As a consequence, the general Majoron cou- 
pling to charged leptons and quarks possesses a scalar and pseudoscalar part. We will 
include these CP-odd effects in our theoretical considerations, although they seem not to 
influence our numerical predictions. 

Armed with the Lagrangians (8)-(10), (14) and (15) it is now straightforward to 
calculate the coupling J" — /i — /2 given by the Feynman graphs shown in Figs. l(a)-l(c). 
The additional CP-violating diagrams of the Majoron coupling to fermions are depicted 
in Figs. 2(d) and 2(e). In our analytical calculations we have neglected terms proportional 
to the small quantities m^/M^ for f = e,u or d. The individual amplitudes contributing 
to the J° — /i — /2 coupling are given by 



-(Aj -|- Aj)(5jj — Cj* )/i(Aj, Xj) + Cij^J AjAj/i(Aj, Xj 



+ -iX..-X,)C:^h{K,X,) 



+ -A^.(A.-A,)C*. /2(A„A,)-/i(A.,A,) 



'-rhh 



{puv - ^) [XiSij - Cij^JXiXj - -{Xi + Xj)C*^ - (Ai + Xj){5i 



^ij)-^'2{^ii^j) + ^ijy^i^ji ~2L2(Aj,Aj) -|- -{Xj — Xi)Is{Xi, Xj 



+ CijL2{Xi, Xj 



(19) 
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(20) 



n-ff 



-i(5^ImCj(cg - setp){Xi - Xj)^J XiXjXf/ Cuv - ^o(Ai, A^) ), 



(21) 



where 



ib^mClAse + cetp)[Xi - Xj)JXiX.jXs'[ Cuv - ^o(Ai, A^) 



:22) 
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A. 
uv 
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A. 
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tan/5 Bi^^Bi^j ui^ 
tan/3 

tan P U/U/ , 
tan /3 uy75U/ . 



Ml 

Mi;' 



q+75) + ^(1-75) 



mil 
Mw' 



Mw 



(23) 
(24) 

(25) 



(26) 



The analytical expressions for the one- loop functions Ji, I2, I3, Lq, Li and L2 are given in 
Appendix A. In Eq. (20) T/ stands for the third component of the weak isospin and takes 
the values: = 1/2, Tl''^ = -1/2. The UV divergences in the amplitudes (19)-(22) 
vanish identically due to the following equalities [jl9|, 



2nG 
1=1 

k=l 



2nr 



J2mnMC:, = 



1=1 

mni^CikCjk 

k=l 



(27) 
(28) 
(29) 
(30) 



Note also that the amplitudes T/-^ , T/^ induce a scalar piece in the — f — f coupling. 



This scalar part, however, is suppressed for astrophysical reasons by a factor of 10 at least 
as compared to the pseudoscalar part of the coupling . 

In order to pin down numerical predictions, we first consider the conservative case of 
a model with one generation or equivalently a three-generation model without interfamily 
mixings. Then, the coupling J° — e — e, gjgg, defined by the relation 



(31) 



takes the simple form 



Sjee 



167r ^M, 



w 



X2 



1 + 



In A 



Ne 



(32) 



The mixings (s^)^ are defined by 



2nG 

E 

i=nG+l 



I Bu I ^ 



(33) 



In this scenario {s^Y = '^']du/^\ii ^^"^ ^ 1, since the heavy Majorana neutrinos 
Ni have to satisfy Eq. (1) and — rni. or rriu^. On the other hand, astrophysical 

constraints arising from helium ignition in red giants or the observational evidence of 
white dwarf cooling rates are given by the bound ^ 



gj,, < (9.-1.4)10 



-13 



(34) 



However, the range 3. 10^^^ < gjgg < 6. 10^^ is excluded from the helium ignition 
argument mentioned above, if the radius of giant core or dwarf is bigger than the mean 
free path of the pseudoscalars that these particles require to freely escape from them |jl|]. 
It is now obvious that the most stringent constraint on tan (3 arises from the heaviest 
family. Thus, for rriD c:^ rrir one obtains that 



3 Jee 



9waw , m 
■ tan p 



32 



MwM^ ' 



(35) 
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yielding because of Eq. (34) 

tan/3 < 0.4. (36) 
Of course, if mo '^nit/k 10 GeV, one finds a much stronger bound, i.e. 

tan/3 < 10"^. (37) 

Note also that such low-energy realizations make unlikely the invisible decay of massive 
Higgses into Majoron pairs f^, |^. 

The afore-mentioned hierarchical scheme, however, is in general not valid if one in- 
troduces intergenerational mixings in the singlet Majoron model. This situation seems 
to be a natural possibilty that can be realized by GUT models, since tud and Mu ma- 
trices may get related in such high-energy scenarios (i.e. m£,{Mx) = Mu{Mx) with Mx 
indicating the grand unification scale). In addition, it has explicitly been demonstrated 



m 



22| , |T9[ that the scale of tum can be C(IOO) GeV without contradicting experimental 
bounds on neutrino masses. For instance, democratic-type mass matrices for the form of 
m£, \2i\ can lead to patterns with such a low scale for niM- Then, the mixings {s^Y can 



be treated as purely phenomenological parameters, since Eq. (33) should now read 

{s^'lf ~ niD^m^D (38) 

mi 

and cannot therefore be related with the light-neutrino mass matrix of Eq. (1). The 
mixing angles {s2Y ^^"^ generally be constrained by a global analysis of a great number 
of low-energy experiments and LEP data P3[. In this scenario one makes the remerkable 



observation that gjgg can severely be supressed for a certain choice of the mass parameters 
\ni and mixings {s^lY- For example, if all heavy neutrino masses uin^ are approximately 
equal and Atv^ ^ 1, then the choice 



sl^f ^ {sr? (39) 
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leads to gjgg = 0. However, even if the Majoron couplings to electrons vanish, the 
corresponding coupling to nucleons TV, gjj\f_\f, is not zero anymore. The reason is that the 
destructive first term in the bracket of Eq. (32) does not exist anymore and such a fine- 
tuning is thus not possible. Since gjj^j^/gjee — i^Nl^^e — 2. 10'^, one may derive useful 
constraints from the consideration of cooling rates of neutron stars due to the energy loss 
mechanism by Majoron emission. In Fig. (3) we present exclusion plots of the parameters 
tan/3 versus tjin for three different values of {s2Y by considering that 



~ 10-^ (40) 



For a discussion of additional uncertaintities on the upper bound of the coupling gj^^ 
that can arise from various reasons like the so-called "Turner's window" |2^ etc., we refer 



the reader to [^]. Ultimately, we must notice that the astrophysical bounds arising from 
the Majoron coupling to two photons, Cj-y-y, should be weaker than that coming from 
SjAfAfy since Cj^^ can only be generated at two-loop electroweak order. 

In the following we will focus our attention on bounds resulting solely from terrestrial 
experiments [|l^, by analyzing lepton- flavor violating decays, i.e. li J°/2 with li 7^/2. 



To the leading order of the heavy neutrino limit one finds from Eqs. (18) and (19) that 

BR(/r ^ J%) ^ ^tan^/? \Blj,B,,r^\'X% (41) 

The experimental information we have for the above lepton-violating decays are the fol- 
lowing upper bounds: 

BR(/i^J°e) < 2.6 lO"*' [15], 
BR(r J°e) < 7.1 10^^ M 



BR(r ^ JV) < 2.3 10"=^ |6|. (42) 



In order to quantitatively estimate the magnitude of the lepton-violating effects that 
could be constrained by the branching ratios stated in (42), we use the upper bound of 
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the quantity 

\BI,,Bi,n\ < is'2r = max((.^)^ {s^^)'). (43) 

The exclusion plots implied by these experiments are presented in Fig. (4) for the three 
different decay channels. For comparison, we have taken the astrophysical bound coming 
from Eq. (40) into account in Fig. (4), from which one easily concludes that experimental 
searches for the decay yU J°e may not be excluded by astrophysical constraints and can 
hence be sensitive to new physics beyond the SM. 

In conclusion, astrophysical considerations may lead to useful constraints on the 
parameters of singlet Majoron models with intergenerational mixings. It has been demon- 
strated that three-generation Majoron models can indeed be constrained if these models 
are assumed to be embedded in GUT scenarios. Possibilities of how to evade from some 
of the astrophysical constraints have also been discussed. For example, gjgg vanishes for 
a specific choice of parameters. Furthermore, terrestrial experiments give independently 
severe restrictions on the lepton-violating mixings and heavy neutrino masses. Aside from 



rather involved i?-parity broken models , this minimal extension of the SM, the singlet 
Majoron model, may also naturally account for possible lepton-flavor violating signals in 
precision experiments. We emphasize again the fact that measurements of the TRIUMF 



collaboration |]T5| for exotic decay modes, like /i J^e, lie in area which may not be 
excluded by astrophysics and have substancial chances to establish new physics beyond 
the SM. Finally, due to the CP-odd interactions that Majorona neutrinos introduce in 
such models (see e.g. Eqs. (9), (14) and (15)), one may be motivated to discuss their phe- 
nomeno logical impact of possible CP- violating effects in the decays of the Higgs particle 
into top, W or Z pairs . 
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A The loop integrals 



We first define the useful functions Bi{Xi, Xj) and B2{Xi, Xj) as: 



Bi{Xi,Xj) = Xi{l-x) + XjX, 
B2{Xi,Xj) = l-y + y[Xi{l-x) + Xjx] , 



(Al) 
(A2) 



where x and y are Feynman parameters. The loop integrals Lq, Li, L2, h, I2 and I3 are 
then given by 



Lo{Xi,Xj) = j dx In Bi{Xi, Xj) 



= -1 + -InAiA, 
2 



+ 



2(Aj — Xj) Aj 



Li(Ai,Aj) = J dxxlaBi(Xi, Xj) 



A. 



4 ■ 2(A,-A,)2 A, 
L2(Aj,Aj) = j dxdy y In B2{Xi,Xj) 



^'"'^ " 2(A.-A,)' 



4(A, - A,; 
1 



A, 



A,; 



1 - A, 1 - A, 



3Aj Aj 



4(1-A,)(1-A,) 



+ 



^AfhiA, 



A?lnA,- 



2(Xi - Xj) [l-Xj l-\ 
r dxdy y 
J Bi{Xi,Xj) 

XiXj ln(Aj/ Xj) + Xj In A^ — Aj In Aj 



72(Aj,A,) = J 



(1-A.)(1-A,-)(A.-A,-) 

dxdy y"^ 



B2{Xi, Xj) 
1 



+ 



1 



+ 



2(1-A,)(1-A,) 2(Aj-A,) 
In Ai 



In 



Aj 



In Aj 



A, (1 - A,)^ 



(1 - A,)2 



^3(Aj,A,.) = / 



dxdy y^{l — 2x) 
B2{Xi, Xj) 
1 

■2(1-Aj)(l-A,) 



In Aj 



In Aj 



1 - Aj 1 - Aj 



2(Aj - A,) 



(A3) 



(A4) 



(A5) 



(A6) 



(A7) 



Aj 



1-A, 
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+ 



A,; 



A, 



ln(A,/A,) 
2(A, - A,)^ 



A, 



+ 



1- A,: 



(A8) 



The integration interval of the variables x and y is [0, 1]. 
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Figure Captions 



Fig. 1: Feynman graphs responsible for the couphng of Majorons to fermions, J° — 

Fig. 2: CP-odd graphs giving rise to a scalar part in the coupling — / — /. 

Fig. 3: Exclusion plots from astrophysical requirements. We have considered the 
values: (s^')^ = 5. IQ-^ (solid line), (s^')^ = 10~^ (dashed line), {s^ f = 
10~^ (dot-dashed line). The area lying above of the curves is excluded by 
the restriction gjxx < 10~^. In addition, we assume that all heavy neutrino 
masses are approximately equal with ttin- 

Fig. 4: Exclusion plots originating from the decays: — > J°e (solid line), r 

J°e (dashed line), r J'^/i (dot-dashed line). For comparison, we have 
considered the astrophysical bound gjj^x < 10~^ (dotted line). The areas 
lying above of the curves are excluded by the afore-mentioned conditions. 
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